Detonation shock dynamics and comparisons with direct numerical simulation

Comparisons between direct numerical simulation (DNS) of detonation and detonation shock dynamics (DSD) is made. The theory of DSD defines the motion of the detonation shock in terms of intrinsic geometry of the shock surface, in particular for condensed phase explosives the shock normal velocity, Dn, the normal acceleration, Ḋn, and the total curvature, κ. In particular, the properties of three intrinsic front evolution laws are studied and compared. These are 1) Constant speed detonation (Huygens’ construction), 2) Curvature dependent speed propagation (Dn − κ relation), and 3) Curvature and speed dependent acceleration (Ḋn −Dn − κ relation). We show that it is possible to measure shock dynamics directly from simulation of the reactive Euler equations and that subsequent numerical solution of the intrinsic partial differential equation for the shock motion (e.g. a Ḋn −Dn − κ relation) reproduces the computed shock motion with high precision. ∗Corresponding Author, Los Alamos National Laboratory, Los Alamos, NM 87545 †University of Illinois , Urbana, IL 61801